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Mode identification for Balloon 090100001 using combined 
multicolour photometry and spectroscopy 
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ABSTRACT 

In this paper, we show that method of mode identification using combined multicolour pho- 
tometry and spectroscopy can be successfully applied to the pulsating subdwarf B star Bal- 
loon 090100001. The method constrains the spherical degree, £.We confirm that the dominant 
mode is radial and we show that for some other modes the method provides values of I con- 
sistent with the observed rotationally split triplet. Moreover, we derive a radius variation of 
1.7 per cent for the dominant mode. The identification opens the possibility for constraining 
the internal structure of the star by means of seismic methods. 

Key words: oscillations - subdwarf - stars: individual: Balloon 090100001. 



1 INTRODUCTION 



About a decade ago, pulsations in hot subdwarfs (sdB) were dis- 
covered i_by_astronomCTsat_the South African Astronomical Obser- 
vatory (Kilk enny et~ai]| 19971 . and subsequent papers in the series). 
At present, there are two classes of pulsating hot subdwarfs (sd- 
BVs) known: V361 Hya and PG 1716 stars. The former pulsate in 
p modes which have short periods, typically a few minutes long, the 
latter, in g modes which have periods an order of magnitude longer. 
Some V361 Hya stars show relatively large amplitudes, reaching or 
even exceeding 50 mmag in V, while for PG 1716 stars only low- 
amplitude modes are observed. The two classes also have slightly 
different effective temperatures, V361 Hya stars are on average hot- 
ter than PG 1716 stars. 

Considering the poorly understood evolutionary history of 
sdB stars, it is very important to obtain information on their in- 
ternal structure whenever possible. The discovery of pulsations in 
sdB stars has opened a new avenue to probe their interiors by means 
of seismic techniques. However, a prerequisite of a successful ap- 
plication of asteroseismology is detection of at least several modes 
which are properly identified in terms of their pulsation geometry. 
The identification can be carried out using different methods that 
utilize photometric or spectroscopic data or a combination of both. 

For sdBV stars, spectroscopic methods have an extra layer 
of difficulty because the currently known objects are fainter than 
V ~ 12 mag, and they all have very short pulsation periods. 
Good time-resolved spectroscopic observations can therefore only 
be obtained using large telescopes. For this reason, there have 
been many attempts to use only (multicolour) photometric data to 
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identify modes in sdBV stars. In particular, methods that compare 
observed and theoretical behaviours of pulsation amplitude with 
wavelength have been used to constrain the spherical degree , I. For 
example, such applic ations have been pr esented by iKoenl d 19981) 
for KPD 2109+440 1, llefferv et all d2004l) for HS 0039+4302 and 
Ijefferv et alj d2005l) for PG 00 14+067. Unfortunately, photometric 
data are limited because the amplitude ratios of modes with dif- 
ferent i sometimes have very similar wave length dependence (see 
iRamachandran. lefferv & To wnsend 2004), and therefore this does 
not lead to an adequate discrimination of I, In the best case, only 
some possibilities could be rejected (see Sect. 5 or Trembl av et al.l 
l2006h . 

Despite the difficulties in obtaining good-quality, 
time-resolved spectra for sdBV stars, spectroscopic time- 
series data have been obtained for some of the bright- 
est objects: for example, PG 1605+072 = V338 Ser 
dO'Toole et aT]|200Cll2002ll2003Ll2005l : IWoolf. Jefferv & Pollaccd 
l2002l: iFalter et al.l 120031) . KPD 2109+44 01 and PB 8 783 
dlefferv & Pollaccol |2OO0h . PG 1627+ 17 dForetalJ l200fj) or 



Balloon 090100001 (hereafter Bal09; iTelting & Ostensenl 120061 
hereafter TO06). In this paper we show that the met hod of 
iDaszvnska-Daszkiewicz. Dziembowski & Pamvatnvkhl d2003l) 
can be successfully applied to the extremely interesting pulsating 
subdwarf Bal09. 



2 BALLOON 090100001 

Bal09 is present ly the brigh t est kn own sdBV star. Its pulsations 
were found by breiro et alj d2004l) from a single night of pho- 
tometric data. The main mode detected in this star has a semi- 
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Figure 1. A schematic Fourier spectrum of the 2004 B-filter data of Bal09. 
Note that semi-amplitudes are given on a logarithmic scale. Modes with the 
largest amplitudes are labelled as in Paper I and Tables 1 and 2. Two com- 
ponents of the quintuplet (not labelled in the bottom panel) were detected 
only in the Fourier spectrum of 2005 data. 



amplitude of about 60 mmag in B and a period of about 6 min, 
quite long for a V361 Hya star. Thus, Bal09 became an obvious 
target for a mor e detailed study. T he follow-up phot ometry car- 
ried o ut in 2004 dBaran et ai]|2005i hereafter Paper I; IOreiro et"ai] 
l2005h revealed several interesting features in its frequency spec- 
trum (Fig. [TJ. First of all, among about 20 detected modes, an 
equidistant triplet was found close to the main mode. If interpreted 
in terms of rotational splitting, it would imply a rotational period of 
about 7 d (Paper I). Even more importantly, several low-frequency g 
modes were found. This made the star one of two presently known 
hybrid V361 Hya/ PG 1716 stars; the other one is HS 0702+6043 
dSchuh et al]2006h . 

In 2005, a large multisite campaign was carried out on Bal09 
(Baran et al., in preparation). The new data allowed the detection 
of a quintuplet (see Fig. []} whose three strongest components had 
already been detected in the 2004 data. In addition, it turned out 
that the amplitudes of some modes change. However, the most sur- 
prising discovery was finding that splittings of multiplets changed 
considerably between 2004 and 2005 (Bara n et all 20071) . As far as 
we are aware, this is the first clear evidence for such a change in 
any pulsating star. This finding poses a problem for the interpre- 
tation of splittings only in terms of rotation. However, we believe 
that rotation is the main factor that causes splittings, while another 
(presently unknown) effect(s) causes its small variations. We there- 
fore refer to the multiplets as rotationally split. 



3 OBSERVATIONS AND ANALYSIS 

The amplitude changes detected in Bal09 and other sdBV stars 
suggest that if one wants to combine photometry and spectroscopy 
for mode identification, the data must be obtained simultaneously. 
For this reason, we decided to use the multiband photometry of Pa- 
per I and spectroscopy of TO06 which were coincidentally acquired 
at the same time, during 2004 August and September. 



The UBVR photometry of Paper I was carried out at two 
sites, Mt. Suhora and Loiano Observatories. About 120 h of data 
were gathered through each filter. The observations were spread 
over about 40 d, resulting in a frequency resolution of better than 
0.5 fiHz. Details of the Fourier analysis of these data can be found 
in the Paper I. As found in that paper, the frequency spectrum of 
Bal09 shows three regions of p modes (around 2.8, 3.8 and 4.7 
mHz), a group of low-frequency (below 0.5 mHz) g modes and 
three regions where combination modes occur, labeled 'comb.' in 
Fig.E 

Because the changes in amplitude appeared to be more pro- 
nounced in the 2005 campaign data than we found in Paper I, we 
decided to reanalyze the 2004 photometry. This time, however, we 
allowed for linear amplitude changes. In addition, the outliers were 
removed with a slightly more stringent criterion than in Paper I. 
This resulted in a lower detection threshold. We also decided to ac- 
cept only those modes which were detected in all four filters. As a 
consequence, the list of modes we present in Tab.[3]is shorter than 
in Paper I. 

As mentioned above, spectroscopic observations of Bal09 
were carried out by TO06, simultaneously with the 2004 photom- 
etry. During seven nights spanning over 38 d, they obtained over 
2500 low-resolution spectra. Importantly the cycle time was only 
43 s long, allowing good sampling of the radial velocity curve. 
Details of the reduction and analysis can be found in TO06. We 
have reanalyzed the radial-velocity data of TO06 in order to obtain 
the fit parameters in a manner consistent with the photometry. As 
the 2004 photometric data suffered from aliasing much less than 
the radial-velocity data of TO06, we always adopted a frequency 
measured from photometry, even if the alias peak was found to be 
higher than that at the "real" frequency. In total, we detected 14 
modes in radial velocity; their parameters are presented in Tab.[T] 
We found the same p modes as TO06: /i, the triplet (/2, fs., and 
fi), fs and two combination frequencies, 2/i and /i + fa. In the 
low-frequency domain, we also detected seven modes above the de- 
tection level. Unfortunately, none of these could be unambiguously 
identified with the modes detected in photometry, and only some 
of these were the same as those reported by TO06. Apparently, the 
rich spectrum of g modes and their very low amplitudes make the 
aliasing a very severe problem during the prewhitening procedure. 

However, three g modes known from photometry, /a, /b and 
fc according to the designations used in Tab. [3] were found by 
TO06 in the variations of the equivalent widths (EW) of the Balmer 
H7-H9 lines. The reason why the same modes were not detected in 
radial velocities can be qualitatively explained if we recall that for 
g modes the horizontal motions dominate the radial motions, while 
for p modes the opposite is true. As all have small amplitudes, it is 
therefore quite possible that different modes are seen in photome- 
try (an d EW changes) and radial velocities. As shown bv lFor et"al] 
(2006), because of the small amplitudes of g modes, it is very dif- 
ficult to detect them in radial-velocity data. For modes observed in 
both photometry and spectroscopy, simultaneous data allow us to 
derive phase lags between light and radial-velocity variations. This 
parameter is related to the non-adiabacity of pulsations. The phase 
lags for modes detected in photometry and spectroscopy are given 
in Tab. [I] The weighted mean of the phase lag for five independent 
p modes is equal to 1.45 ± 0.03 rad or 83 ± 2°. For the dominant 
mode, /1, it is very similar, 1.43 ± 0.02 rad. Because this mode is 
radial (see Sect. 5), the phase of maximum brightness corresponds 
roughly to the phase of minimum radius. In other words, we can 
conclude that in V361 Hya stars the brightness variations are dom- 
inated by the temperature effects. Adopting projection factor equal 



Table 1. The results of re-analysis of 2004 UBVR data of Bal09. Semi-amplitudes, A, and phases, <j>, are given for the epoch HJD = 2 453 250. 
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Table 1. Results of the sine-curve fit to the low-resolution spectroscopic 
data of TO06. As for the photometric data (Tab.[5), the phases are given for 
the epoch HJD = 2 453 250. A<f> = 0rv — </>phot denotes difference between 
radial-velocity and photometric (averaged over all four filters) phases. 
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to 1.4 as typical for an sdB star (see Montafies Rodrigu ez & Jeffervl 
l200ll) and a radius of 0.25 Rq, we estimate the full range of radius 
variations for the dominant mode in Bal09 to be about 1.7 per cent. 



4 THE METHOD OF MODE IDENTIFICATION 

Having derived frequencies, amplitudes and phases 
we can now apply the metho d (described in detail by 
iDaszv hska-Daszkiewi cz et al to discriminate I parameters 

for the modes detected in our photometry and spectroscopy of 
Bal09. 

In the simplest form , the method of 
IDaszvhska-Daszkiewicz et al. I J2003h utilizes a set of linear 
equations, of the form 



(1) 



to derive two complex parameters, e and /. The coefficients A\ are 
taken from observations; these are amplitudes, a\, and phases, <j>\, 
written in a complex form, A\ = a\ exp(i(f>\), where A denotes a 
given photometric band. However, D\j and £x,e depend on the the 
flux and limb darkening derivatives over effective temperature and 
surface gravity, and have to be calculated using model stellar at- 
mospheres. In principle, observations in two photometric bands are 
sufficient to solve the set of equation (1) once the spherical degree 
I is assumed. However, if observations in more than two photomet- 
ric bands are available, the set becomes overconditioned and can be 
solved by means of the least-squares method. Usually, however, we 
are not interested in the values of e and /, but in discriminating I. 
In order to achieve this, the equations are solved for different val- 
ues of £ and then a goodness-of-fit statistic like reduced \ 2 is used 
to decide on discrimination. 

Unfortunately, as was shown 
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see also Sect. 5), using merely photometric data rarely leads to 
unambiguous discrimination of £. It is much better to combine 
photometric and spectroscopic data adding another equation of set, 
that is 



Tx,ie = Mi 



(2) 



In Eequation (2), M\,\ represents the first moment of spectral line 
variations, (i.e. radial velocities), while T\,t is again calculated 
from the model atmospheres. The solution of a set of equations 
(1) and (2) leads in this case to a much better discrimination of i 
(see Section 5). 

The model-dependent coefficients in equations (1) & (2) and, 
consequently, the discrimination of £ depend on the model parame- 
ters, in particular effective temperature, surface gravity and metal- 
licity. For B al09, the effective te mperature and surface gravity were 
derived by lOreiro et alj J2004h. T ff = 29 450 ± 500 K, log g = 
5.33 ± 0.10 dex and iTeltinget alJ feOOd) . T cff = 28 700 K, log g = 
5.39 dex. We therefore calculated flux distributions in the range be- 
tween 300 and 1000 nm for a grid of models with T e g ranging from 
28 100 and 31 000 K and logs from 5.2 to 5.5 dex. To do this, a 
grid of fully metal line-blanketed local thermodynamic equilibrium 
(LTE) s tellar atmosphe re models was calculated with the ATLAS9 
code of lKurucj A1992). Flux distributions were then calculated us- 
ing Michael Lemke's version of the LINFOR program (originally 
developed by Holweger, Steffen, and Ste enbock at Kiel University ), 
similar to the synthetic spectra used bv lO'Toole&He ber (2006). 
The Kurucz line list was used as a source of oscillator strengths 
and damping constants for all metal lines. Only lines that have 
been observed experimentally were included. The resulting flux 
distributions were multiplied by transmission functions of the filters 
in order to obtain integrated fluxes in the four photometric bands, 
UBVR, that we used. 

For the same grid of models, we also calculated the specific in- 
tensities as a function of /j, = cos 8, where 8 was the angle between 
the line of sight and the normal to the stellar surface in a given point 
of the stellar disc. Using these intensities, the limb darkening law 
was derived, fitting the coefficients c and d of a function 

JO*) = /(l)[l-c(l-/i)-d(l-VJi)] O) 

dPfaz-Cordoves & Gimenezlll992l) which was found to reproduce 
best the calculated changes of I(fi). This was carried out for differ- 
ent combinations of input parameters of the model, thus allowing 
the calculation of neccessary derivatives. 

The flux distribution calculated from a model atmosphere de- 
pends also on helium abundance, metallic ity and mi c roturb ulent 
velocity. Following the determination of lOreiro et alj d2004l) . he- 
lium abundance was assumed to be log[n(He)/n(H)] = —2.54. In 
addition, we assumed solar metallicity and microturbulent veloc- 
ity £ = kms -1 . The latter two parameters have not been deter- 
mined yet for Bal09, b ut are typical values for an sdB star (e.g. 
lO'Toole & Heber 20(j3). 



5 DISCRIMINATION OF THE I PARAMETER 

As a first step, we used the method presented in Section 4 in its 
simplest form (i.e. taking into account only photometric amplitudes 
and phases in four bands; see Tab. 0. The results for six modes 
detected both in photometry and radial velocities are presented in 
Tab.|2]and in the left panels of Fig.[2] We limited our calculations to 
£ ^ 4 because, as a result of averaging over the stellar disc, modes 
with larger I are not expected to be visible in photometry. 

As can be seen from Fig. [2] and Tab. [2] there is only little dis- 
crimination of I for /1, /2 and f$ and none for ft and /s. For /1, £ 
= or 2 are possible, while for fa and fa, I > 2 are excluded. 

The poor discrimination of I when using photometric data can 
be understood when considering t he amplitude ratios ve r sus w ave- 
length calculations carried out by Ramachandra n et al.1 d2004l) . In 



Table 2. The results of the application of the method of 
Daszyhska-Daszkiewicz et al ] J2003I) for Bal09 2004 photometric data: 
reduced x 2 values for five modes with the largest amplitudes as a function 
of the spherical degree I, 



I 


/i 


h 


h 


h 


fa 





3.56 


1.10 


0.81 


0.48 


0.26 


1 


13.22 


1.03 


0.89 


0.84 


0.28 


2 


7.48 


1.09 


0.89 


0.57 


0.18 


3 


42.12 


10.67 


5.21 


0.46 


0.43 


4 


28.39 


5.07 


2.30 


0.34 


0.46 
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Figure 2. Reduced \' 2 as a function of the spherical degree £ for 2004 data 
of Bal09 and five modes detected both in photometry and spectroscopy. The 
modes are labelled with the encircled designation in the middle (see Tables 
[5]and[T}. Left panels show the results when only photometric data were 
used in calculations, and right panels when photometric and spectroscopic 
data were combined. The i value(s) that we regard as possible for a given 
mode are shown with large circles. 

addition, three of the four photometric bands we used are located in 
the Paschen continuum. The resulting equations for different bands 
might therefore be highly correlated, which will not lead to a good 
discrimination of £. 

It also has to be pointed out that the reduced \ 2 statistics we 
use to discriminate £, does not have, in our case, its proper char- 
acteristics, as the set of equation (1) is not well overconditioned; it 
has only four degrees of freedom. For this reason, we did not set a 
given threshold in to indicate possible values of £. Instead, we 
do this arbitrarily, comparing the values of \ 2 f° r different values 
of £; those with the smallest % 2 are regarded as possible (see Fig.O. 

In the next step, we added the radial velocity data as an addi- 
tional constraint via equation (2) and repeated the procedure. The 
results are shown in the right panels of Fig. [2] and in Tab. [3] As 
equation (2) contains only e as a free parameter, it is expected that 
/ and e will be much better constrained than previously and, con- 
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Table 3. The same as in Tab. [2] but for the combined photometric and spec- 
troscopic data. 



£ 


/l 


h 


h 


h 


h 





3.88 


4.44 


3.03 


0.96 


0.22 


1 


9.03 


0.80 


0.94 


0.65 


0.35 


2 


47.73 


1.99 


1.06 


0.51 


0.55 


3 


218.07 


77.42 


28.11 


4.63 


3.56 


4 


1526.90 


101.09 


35.45 


6.03 


6.03 



Table 4. The possib le values of £ for five s trongest modes in Bal09 derived 
in this paper and by Tremblav et al. 1 2006). 



Mode 


This paper 


This paper 


Tremblav et al. t2006) 




(photometry) 


(photometry + 








spectroscopy) 




/i 


0,2 





0, 1 


h 


0, 1,2 


1 


1,2 


h 


0,1,2 


1,2 


1,2 


h 


0, 1,2,3,4 


0,1,2 


1,2,4 


h 


0, 1,2,3,4 


0, 1,2 


0,1,2,3 



sequently, £ will be much better discriminated. As can be seen from 
Fig. [2] this is indeed the case. Adding spectroscopy improves the 
situation in the sense that observational equations become less cor- 
related, which allows reliable estimation of e and /, but not as far as 
the x 2 statistics is concerned. For combined data, the problem has 
two more (i.e. six) degrees of freedom, still too small a number to 
say the problem is well defined. Because of that, we would like to 
avoid giving a level of significance, as it would be irrelevant. Nev- 
ertheless, relative numbers tell us which values of £ are acceptable, 
although, as we have pointed out, our choice is arbitrary. 

The best and least ambiguous discrimination was obtained for 
the dominant mode, which appears to be radial (£ = 0). For the 
strongest component of the triplet, / 2 , £ = 1 is definitely the best 
solution, exactly as expected. For the two remaining triplet com- 
ponents, fs and /4, £ is not unambiguously identified, but £ = 1 
is one of the possibilities. Therefore, we conclude that, as already 
suggested in Paper I, the main mode in Bal09 is radial, while the 
triplet represents a rotationally split 1=1 mode. 

Table [4] summarizes the results of our discrimination. In 
addition, we compar e these results with those presented by 
iTremblav et alj ([2006), who used only photometric data to iden- 
tify £ for the nine modes published in Paper I. They did not obtain 
a unique identification for any of these modes. For the best cases, 
they present two possible values of £, similar to the results we ob- 
tained using only the photometric data. 

We also checked how the discrimination we obtain depends on 
the parameters of the model. The result for adopting five different 
effective temperatures from the range between 28 500 and 30 500 K 
is shown in Fig. [3] The change with T e g does not greatly affect the 
result of discrimination, although lower effective temperatures £ = 
1 for /i and £ = 2 for fa become alternatives to those given in 
Tab. [4] A negligible effect was obtained while changing the mass 
between 0.47 and 0.53 Mq and radius between 0.23 and 0.26 Rq. 
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Figure 3. Discrimination of the i parameter for five different values of T c g, 
combined photometry and spectroscopy and two strongest modes, /i (top) 
and /2 (bottom). Going from the shortest to the longest dashes results are 
shown for effective temperatures between 28 500 and 30500 K in 500-K 
steps. 



6 SUMMARY AND CONCLUSIONS 

Using photometric and spectroscopic data for Bal09, we have 
shown that the method of mode identification developed by 
iDaszvriska-D aszkiewi cz et al. I J2003I) can be successfully applied 
to Bal09. We believe that, if simultaneous multicolour photometry 
and spectroscopy are carried out, the method will also allow mode 
identification in other sdBV stars. In Bal09, the spherical degree 
£ was unambiguously identified for the two modes with the largest 
amplitudes. For several other modes, we limited the number of pos- 
sible values of I to two or three. If we accept that the two observed 
multiplets are caused by rotational splitting, then the number of 
identified modes increases to nine: the triplet components have t = 
1, the quintuplet components, 1 = 2. 

It is clear from our considerations that photometry itself is not 
sufficient to identify mode satisfactorily for sdBV stars, and that 
adding spectroscopic observations greatly improves the situation. 
For Bal09, it is certainly worth carrying out another photometric 
and spectroscopic campaign to identify modes with even lower am- 
plitudes. However, our identification is probably sufficient to go a 
step further and attempt asteroseismology of this star. However, for 
this purpose, an appropriate set of models is req uired. A first at- 
tempt was presented in Paper I using the models of Charpinet et al. 
(2002). While we had no modes identified at that time, the assump- 
tion that the main mode is radial now appears to be correct. With 
reliable identification of modes from the 2.8 mHz group, a good 
seismic model that fits the identified modes might be used to match 
(and identify) remaining modes. In order to carry out a thorough as- 
teroseismological analysis and potentially learn new physics, it is 
important that more modes are identified independent of pulsation 
models in sdBV stars. We have provided a crucial piece of the puz- 
zle, and the challenge is now for pulsation theorists to find models 
that match both the pulsation frequencies and mode identifications 
derived here. 
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